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Topological order characterizes those phases of matter that defy a description in terms of symmetry and 
cannot be distinguished in terms local order parameters. This type of order plays a key role in the theory of 
the fractional quantum Hall effect, as well as in topological quantum information processing. Here we show 
that a system of n spins forming a lattice on a Riemann surface can undergo a second order quantum phase 
transition between a spin-polarized phase and a string-net condensed phase. This is an example of a quantum 
phase transition between magnetic and topological order. We furthermore show how to prepare the topologically 
ordered phase through adiabatic evolution in a time that is upper bounded by O(y'n). This provides a physically 
plausible method for constructing and initializing a topological quantum memory. 
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The notion of topological order yj] can explain phase tran- 
sitions that depart from the standard Landau theory of symme- 
try breaking, long range correlations and local order parame- 
ters. It plays a key role in condensed matter theory of strongly 
correlated electrons; e.g., it provides a means to understand 
the different phases arising in the fractional quantum Hall ef- 
fect |2[]. These phases have exactly the same symmetries. No 
local order parameter can distinguish them. The internal or- 
der that characterizes these phases is topological, and can be 
characterized in terms of a non-local order parameter, e.g., the 
expectation value of string operators (operators formed by ten- 
sor products of local operators taken along a siring). Recently, 
it was shown that topological order can also be characterized 
by topological entropy yfl. 

Another arena in which topological order has found pro- 
found applications is quantum computation, in particular in 
the context of systems exhibiting natural fault tolerance 
Central to this application is the ability to prepare certain 
topologically ordered states, which are the ground states of a 
Hamiltonian describing a spin system occupying lattice links 
on a Riemann surface of genus g (e.g., a torus, for which 
g = 1). In a topologically ordered phase the ground state 
degeneracy depends on g. For instance, in the toric code [14] 
the ground state manifold C is 2 2fl -fold degenerate, and this 
allows one to encode 2g qubits. Because any two orthogonal 
states in C are coupled only by loop operators with loops that 
wind around the holes of the Riemann surface, this encoding 
is robust against any local perturbation. Such a system consti- 
tutes a very robust memory register 14J] . 

Here we address the problem of preparation of topologi- 
cally ordered states via adiabatic evolution. Our motivation 
for studying this problem is at least threefold. First, while 
quantum phase transitions (QPTs) between deconfined and 
confined phases have been studied extensively in quantum 
chromodynamics [5], this has not been the case in condensed 
matter physics. Here we study a second order QPT between a 
magnetically ordered phase (deconfined, topologically disor- 
dered) and a string-net condensed phase [8] (confined, topo- 
logically ordered). Second, we are motivated by the afore- 
mentioned need for topologically ordered states as the input 



to topological quantum computers. In Ref. fioll it was shown 
how to prepare a ground state of the toric code via 0(y/n) 
repeated syndrome measurements, where n ~ L 2 is the num- 
ber of spins on a lattice of linear dimension L. Such ground 
state preparation is an essential step in constructing a topo- 
logically fault tolerant quantum memory. However, prepara- 
tion via syndrome measurements has certain drawbacks, most 
notably that the measurements are assumed to be as fast as 
logic gates and must be fast compared to the decoherence 
timescale, assumptions which are likely to be hard to satisfy 
in practice. A system in which the phase transition predicted 
here can potentially be realized experimentally, as well as 
used for topological quantum memory, is a Josephson junc- 
tion array jl 111 . Third, we believe that the methods presented 
here will also find applications in the field of adiabatic quan- 
tum computing [12], for the preparation time in our adiabatic 
method is 0{^/n) (as in the syndrome measurement method 
ifioll ). which is optimal in the sense that it saturates the Lieb- 
Robinson bound II 1311 . 

The toric code model. — Consider an (irregular) lattice on 
a Riemann surface of genus g. At every link of the lattice we 
place a spin 1/2. If n is the number of links, the Hilbert space 
of such a system is U = Hf n , where Hi = Span{|0), |1)} is 
the Hilbert space of a single spin. In the usual computational 
basis we define a reference basis vector 1 0) = 1 0) i (g> . . . <g> 1 0) „, 
i.e., all spins up. This string-free state is the vacuum state for 
strings. We define the Abelian group X as the group of all 
spin flips on H. Its elements can be represented in terms of 
the standard Pauli matrices as x a = {cr'f) ai ® ... ® (c^) Q ", 
where ay S {0, 1} and a = {<x/}™ =1 . Every x a S X squares 
to identity I, dim Tt = \X\ = 2™, and any computational basis 
vector can be written as |a) = x a \0). 

Now we associate a geometric interpretation with the ele- 
ments of X, as in Fig.Q] To every string 7 connecting any two 
vertices of the lattice we can associate the string operator x 1 
operating with a x on all the spins covered by 7. The product 
of two string operators is a string-net operator: x 1 = x 11 x l2 , 
and we can view X as the group of string-net operators on 
the lattice. A particularly interesting subgroup of X is formed 
by the products of all closed strings; we denote this closed 
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FIG. 1: (color online) The square lattice on a torus. Closed a>strings 
(in red) commute with the star operator A s because they have an 
even number of links in common. The elementary closed 2>string is 
the plaquette B p . Dashed z-strings connect plaquettes p. The two 
incontractible s-strings are denoted ti,2. The dual string variable 
(dashed, blue) connects the reference lines t\ , with the vertices p 
in the dual lattice. Thus the dual operator anti-commutes 
with the plaquette variable fi x (p) = B v . 



string-nets group jl] by X. If the Riemann surface has genus 
g, we can also draw 2g incontractible strings around the holes 
in the surface: T — (t\, t^g); see Fig. [TJ The group of 
contractible closed strings is denoted by B, and is generated 
by the elementary closed strings: B p — Yijedp a j' where 3 
labels all the spins located on the boundary dp of a plaquette 
p. A star operator is defined as acting with a z on all the spins 
coming out of the vertex s: A s = Yijes a j- Then X splits 
into 2g topological sectors labeled by the tf. X = T ■ B 11411 . 

Topological order in the toric code model. — Consider the 
Hamiltonian H(g,U) = -gJ2 P B P ~ ^E s ^- For g = 
U = 1 this is the Hamiltonian introduced by Kitaev fl, and 



we will refer to it as the Kitaev Hamiltonian. It is an example 
of a zVspin liquid, a model that features string condensation 
in the ground state. The operators A s , B p all commute, so the 
model is exactly solvable, and the ground state manifold is 
given by C = -j» G H | A s \<j>) = B p \<j>) = \<p)}. The vector 
in the trivial topological sector of C is explicitly given by 11411 

i</> > = iBr*$>i°>- (i) 

x£B 

On a torus the ground state manifold C is 4-fold degener- 
ate. A vector in C can be written as the superposition of the 
four orthogonal ground states in different topological sectors 

+ 2j. On an arbitrary 



t\ti\<f>o),i : j € {0,1}, k = 
lattice (regular or irregular) on a generic Riemann surface of 
genus g there are 2g such operators ti and thus the degener- 
acy of the ground state manifold is, in general 2 2fl . This is a 
manifestation of topological order Jl]]. 

A quantum phase transition. — Consider the Hamiltonian 
obtained by applying a magnetic field in the z -direction to all 
the spins: 



£,nl. 



(2) 



The zero-energy ground state of this Hamiltonian is the mag- 
netically ordered, spin-polarized string-vacuum state |0). In 
the presence of a string state x 7 10), with x 1 — Yljej a j ^ 
I, pays an excitation energy that depends only on the string 
length Z 7 , namely (0|a: 7 i?^j: 7 |0) = 2^ 7 . Thus is a tension 
term. If we add the star term Hy = — U J2 S A s to H^, the 
ground state is still |0) [its energy is E g (U) = —Un s ], but the 
spectrum is quite different. Now there is a drastic difference 
between open and closed x-type strings. Closed strings com- 
mute with all the A s and hence only pay the tension, whereas 
open strings also pay an energy of 2U for every spin that is 
flipped due to anti-commutation (—U(0\ajA s aj\0) = +U). 
So every open string pays a total energy of AU. In the U ^> £ 
limit, open strings are energetically forbidden. On the other 
hand, the plaquette operators B p deform ^-strings: the pla- 
quette term H g = —g ^2 p B p acts as a kinetic energy for 
the (closed) strings and induces them to fluctuate. In light 
of these considerations the total model H = + H(g, U) 
has two different phases. The first phase is the spin polar- 
ized phase described above, when U ^ £ ^S> g ~ 0. Here 
the string fluctuations are energetically suppressed by the ten- 
sion, so the ground state is |0). The other phase arises when 
U ^> g £ ~ 0. Now the tension is too weak to pre- 
vent the closed loops from fluctuating strongly. Nevertheless, 
open x-type strings cost the large energy AU ', so are forbidden 
in the ground state. The ground state consists of the super- 
position with equal amplitude of all possible closed strings, 
C. This is the string-net condensed phase. As the tension 
decreases, and the fluctuations increase, the ground state be- 
comes a superposition of an increasingly larger number of 
closed strings. Open strings cannot be present in the ground 
state because their energy is too large. In the thermodynamic 
limit, for some critical value of the ratio between tension and 
fluctuations £/g, the gap between the ground state and the 
first excited state closes and the system undergoes a second 
order QPT. Notice that the gap between the ground and first 
excited state is given only by the interplay between tension 
and fluctuations, and is unaffected by Hij. The second phase 
we have described cannot be characterized by a local order 
parameter, such as magnetization: we have phases without 
symmetry breaking. This is an example of topological order. 
The low energy sector has energy much smaller than U and it 
thus comprises closed strings: £i ow = span{x|0), x S X}. 
This subspace also splits into four sectors labeled by the ele- 
ments of T: = span^i^l )^ G = 0,1}. Thus 
dimAow = dimCH)/(2"/ 2+1 ). Within each sector, low level 
excitations have an energy E\ ow < g and become gapless only 
at the critical point in the thermodynamic limit. 

Adiabatic evolution. — We now show how to prepare the 
topologically ordered state |</>o) through adiabatic evolution. 
The idea is to adiabatically interpolate between an initial 
Hamiltonian H(Q) whose ground state is easily preparable, 
and a final Hamiltonian H{1) whose ground state is the de- 
sired one. The interpolation between two such Hamiltonians 
has been used as a paradigm for adiabatic quantum computa- 
tion 11211 . This process must be accomplished such that the 
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error S between the actual final state and the desired (ideal 
adiabatic) ground state of H(l), is as small as possible. The 
problem is that real and virtual excitations can mix excited 
states with the desired final state, thus lowering the fidelity 
111 511 . Rigorous criteria are known and used below, 

which improve on the first-order approximation typically en- 
countered in quantum mechanics textbooks. For smooth inter- 
polations and one relevant excited state it has been proven that 
an arbitrarily small error 6 is obtained when the evolution time 
T obeys T = 0(1/ AE m i n ), where A£J m ; n is the minimum 
energy gap encountered along the adiabatic evolution [18]. 
Consider the following time-dependent Hamiltonian: 



H(t) =H v + [1- /(t)] Hz + f( T )H g 



(3) 



Where U > 5, £ - 1 and / : r e [0,1] 1— »• [0,1] such 
that /(0) = and /(l) = 1, r = t/T is a dimensionless 
time. At t = the Hamiltonian is H(Q) = Hjj + H^ whose 
ground state is the spin polarized phase |0). At r = 1 the 
Hamiltonian is the Kitaev Hamiltonian H (1) = Hu + H g 
whose ground state is in C, the string-net condensed phase. 
The tension strength is given by £(1 — /(r)), while the fluc- 
tuations have strength gf(r). Adiabatic evolution from the 
initial state |0) will cause the fluctuations to prepare the de- 
sired ground state \4>o). Notice that because of the interplay 
between the tension term H^ and the large term Hu, at ev- 
ery r the ground state is in C^ w . As the tension decreases, 
the gap between the ground state and states in the other sec- 
tors C[ 3 ow (i + j > 0) diminishes. Beyond the critical point, 
the gap scales down exponentially with the number of spins. 
Nevertheless, this does not cause transitions: irrespective of 
how small the gap is, the Hamiltonian does not couple be- 
tween different sectors. Indeed, the adiabatic theorem states 
that the probability of a transition between the instantaneous 
eigenstate |V'j( T )) the system occupies and another eigenstate 
\ipj(r)) (with respective energies Ei(r),Ej(r)) is given by 
Pij < \{ipi{T)m 3 iT))/{Ei{T) - £,-(r)) 2 | 2 . In our case, 
(^i(r)l-H'l^-(r)) = whenever \ipi(r)), belong to 

two different sectors 1I20I1 . Thus transitions between different 
£ low sectors are completely forbidden, and the adiabatic evo- 
lution keeps the ground state in the sector it starts in. The 
situation is different in the presence of an arbitrary fc-local 
perturbation V: time evolution can then, in principle, couple 
the different sectors, and this requires a careful analysis. Be- 
fore topological order is established different sectors remain 
gapped because of the tension, and the adiabatic criterion ap- 
plies. Once in the topologically ordered phase, the exponen- 
tially small gap requires us to apply time dependent perturba- 
tion theory to compute the probability of tunnelling between 
different sectors. It turns out that only the i/fcth order in the 
Dyson series is non-vanishing, and the result is that as long 
V < £, these transitions are suppressed exponentially in L/k 
at arbitrary times 12011 . Namely, the probability of transition 
between two different sectors a and b is 



This result is confirmed by numerical analysis l2lll : topolog- 
ical order protects the adiabatic evolution from tunneling to 
other sectors. Note that, in contrast, preparing the ground 
state by cooling would result in an arbitrary superposition in 
the ground state manifold. Therefore we are concerned with 
the gap within C^ w . We next show that this gap scales as 
AE(n]T c ) ~ n -1 / 2 at the critical point r c . To do this, we 
must understand the lattice gauge theory that represents the 
low energy sector of the model. 

Lattice gauge theory and mapping to the Ising model. — 
We briefly review the lattice gauge theory invented by Wegner 
01911 . showing that the low energy sector (E <C U) of Kitaev's 
model maps onto this theory. The lattice gauge Hamiltonian 
is given by 



s,k 



k)-X; 



(5) 



A local gauge transformation consists in flipping all the spins 
originating from s, so it can be defined as: A s = Yl 



ies a j- 

It is simple to check that H gaugc is invariant under this local 
gauge transformation. Elitzur's theorem [23] implies that the 
local symmetry cannot be spontaneously broken. In a gauge 
theory the Hilbert space is restricted to the gauge-invariant 
states. Thus the Hilbert space for this model is given by 
H g augo = {M G U I A s \ip) = \ip)} C H, the Hilbert space 
spanned by all possible spin configurations. The Hamiltonian 
(0 clearly resembles our H(t), apart from the absence of the 
star term Hjj. Consequently Ti gaugc is just the low energy 
sector E <C U of H (t). In the limit U — > 00 the Hamiltonian 
H(t) maps to # ga u g e with A2/A1 = g/(r)/£(l - /(t)). For 
a critical value of X1/X2 ~ 0.43 H2lL 12211 . the system under- 
goes a second order QPT 

The lattice gauge theory is dual to the 2 + 1-dimensional 
Ising model [12411 . To show this we define a duality mapping, 
as follows. To every plaquette of the original lattice we as- 
sociate a vertex p on the new lattice, and we define the first 
dual variable p, x (p) = B p . In order to obtain the correct 
mapping, the second dual variable must realize the Pauli al- 
gebra with n x (p). Consider now the string 7_ > («(p) from 
the reference line t^ti) to the vertex p on the new lattice, 
see Fig. [T] To this string we associate the operator that con- 
sists of a z on every link intersected by 7^(|)(p) and denote 
it by H z ^^{p). This is the second dual variable. These vari- 
ables realize the Pauli algebra: (n x (p)) 2 = (l^ z (p)) 2 = 1 
and {fi x (p'), n z (p)} + = Sp'p. In order to write H gaugc in 
terms of the dual variables, we note that ^(p)fi^(p — y) = 
<j z (s,x) and n z _^{p) ^{p — x) = cr z (s,y). Then the map- 
ping (a x ,(j z ) 1— > (/i x ,/i z ) yields 



a. 



H 



Ising 



-A 



2^M J (P)-Al 
P 



fi z (p)n z (p+i) 7 (6) 



WWr) < \(V/0 L/k L\ 



Vt 



(4) 



where /1 2 is of the — *,t type if i = x,y respectively. We 
recognize this Hamiltonian as the quantum Hamiltonian for 
the 2 + 1-dimensional Ising model [24]. This model has two 
phases with a well understood second order QPT separating 
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them. The gap between the ground state and the first excited 
state of this model are known to scale at the critical point 
as A(L) ~ L^ 1 where L — n 1 / 2 is the size of the system 
II25I1 . The first excited state is non-degenerate. Returning to 
our H(t), and recalling that its low energy sector E <C U 
corresponds to the spectrum of the gauge theory Hamiltonian 
(0, it now follows that also the gap of H(t) to the first ex- 
cited state scales as A(n) ~ nT 1 ' 2 near the critical point. 
Hence we know that for a smooth interpolation the adiabatic 
time T = 0(1/A£ m i n ) scales as n 1 / 2 , with the final state 
arbitrarily close to |0o) 11811 . 

Error estimates. — The rest of the spectrum consists of two 
large bands with gaps of order g, U. We now wish to estimate 
how well the actual final state \tp(t = T)) [the solution of 
the time-dependent Schrodinger equation with H(t)] approx- 
imates the desired adiabatic state |0o) at t = T [the instan- 
taneous eigenstate of H(l)], given that there is "leakage" to 
the rest of the spectrum. To this end we use the exponential 
error estimate 11611 : for a continuously differentiable Hamil- 
tonian H(t), S = ||V>(T) - 0o COU < C(n) exp [— A(n)T], 
where A(n) is the relevant gap and C(n) ~ ||0||/A(n), with 
4>(t) the instantaneous eigenstate. Since g <C U, transitions 
into the lower band dominate. In this case we can show that 
|| (/>|| ~ 0(n) and hence that already an adiabatic time as short 
as T = (p/g) logn gives an arbitrarily small error S < n 1-p 
(forp > 1) [20]. Therefore the true adiabatic timescale is set 
in our case by the previously derived T = (^(n 1 ' 2 ), which 
completely suppresses errors due to leakage to other gapped 
excited states. We have thus shown that the topologically or- 
dered state |0q) can be prepared via adiabatic evolution in a 
time that scales as n 1 / 2 . This is consistent with the strategy 
devised in [ 10], where a toric code (the ground state |0o)) is 
prepared via a number of syndrome measurements of order 
L = n 1 / 2 . This scaling has recently been shown to be op- 
timal using the Lieb-Robinson bound 11311 . thus proving that 
the adiabatic scaling found here is optimal too. 

Outlook. — In this work we have shown how to prepare a 
topologically ordered state, the ground state of Kitaev's toric 
code model, using quantum adiabatic evolution. The adia- 
batic evolution time scales as (^(n 1 / 2 ) where n is the number 
of spins in the system. This method gives a measurement-free 
physical technique to prepare a topologically ordered state or 
initialize a topological quantum memory for topological quan- 
tum computation H@,0]. 

We can extend Kitaev's model to 3D. In this case, the star 
operators A s generate surfaces called membranes. It has been 
shown that this model too has topological order |9[]. Its low- 
energy sector maps onto the 3D quantum lattice gauge the- 
ory, whose free-energy expansion shows a first order QPT 
between a confined and deconfined (membrane-condensed) 
phase [26]. We therefore conjecture that adiabatic evolution 
from a magnetically ordered phase to a membrane-condensed 
one |@] will involve a first order QPT. 

The robustness of the ground state of the Kitaev model to 
thermal excitations Q4f| endows our preparation method with 
topological protection against decoherence for t > T; how 



to extend this robustness to all times t is another interesting 
open question, the study of which may benefit from recent 
ideas merging quantum error correcting codes and dynamical 
decoupling with adiabatic quantum computation 12711 . An ex- 
ample of the robustness of the topological phase studied here, 
when the spin system interacts with an ohmic bath, has re- 
cently been provided in Ref. [28]. 

Finally, we note that all problems in the computational class 
"Statistical Zero Knowledge" (SZK - a class which contains 
many natural problems for quantum computers), can be re- 
duced to adiabatic quantum state generation 12911 . and prepara- 
tion of topological order may yield new topological problems 
in SZK. 
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